Future CMB cosmological constraints in a dark coupled universe 
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Cosmic Microwave Background satellite missions as the on-going Planck experiment are expected 
to provide the strongest constraints on a wide set of cosmological parameters. Those constraints, 
however, could be weakened when the assumption of a cosmological constant as the dark energy 
component is removed. Here we show that it will indeed be the case when there exists a coupling 
among the dark energy and the dark matter fluids. In particular, the expected errors on key 
parameters as the cold dark matter density and the angular diameter distance at decoupling are 
significantly larger when a dark coupling is introduced. We show that it will be the case also for 
future satellite missions as EPIC, unless CMB lensing extraction is performed. 

PACS numbers: 



I. INTRODUCTION 

Current cosmological measurements point to a flat 
universe whose mass-energy includes 5% ordinary matter 
and 22% non-baryonic dark matter, but is dominated by 
the dark energy component, identified as the engine for 
the accelerated expansion [IH3I- The most economical 
description of the cosmological measurements attributes 
the dark energy to a Cosmological Constant (CC) in 
Einstein's equations, representing an invariable vacuum 
energy density. The equation of state of the dark energy 
component w in the CC case is constant and w = — 1. 
However, from the quantum field approach, the bare 
prediction for the current vacuum energy density is 
~ 120 orders of magnitude larger than the observed 
value. This situation is the so-called CC problem. In 
addition, there is no proposal which explains naturally 
why the matter and the vacuum energy densities give 
similar contributions to the universe's energy budget at 
this moment in the cosmic history. This is the so-called 
why now? problem, and a possible way to alleviate 
it is to assume a time varying, dynamical fluid. The 
quintessence option consists on a cosmic scalar field 
<f> (called quintessence itself) which changes with time 
and varies across space, and it is slowly approaching its 
ground state. Also, the quintessence equation of state is 
generally not constant through cosmic time 
principle, the quintessence field may couple to the other 
fields. In practice, observations strongly constrain the 
couplings to ordinary matter (l4| . However, interactions 
within the dark sectors, i.e. between dark matter 
and dark energy, are still allowed. This could change 
significantly the universe and the density perturbation 
evolution, the latter being seeds for structure formation. 
For models equivalent to the one studied here, see e.g. 
Refs. EErJU. 



affect the cosmological constraints expected from future 
CMB experiments. The Planck satellite mission, for ex- 
ample, is expected to provide high quality constraints on 
several key parameters (see e.g. [23|]). However, those 
forecasts are usually performed under the assumption 
that the dark energy component is either a cosmologi- 
cal constant or a fluid with constant, redshift indepen- 
dent equation of state w = P/p. It is therefore timely 
to investigate if the assumption of a more elaborate dark 
energy component with a coupling with the dark mat- 
ter could have an impact on these constraints. Here we 
indeed focus on the future CMB data constraints on in- 
teracting dark matter-dark energy models, exploiting, in 
particular, the gravitational CMB lensing signal. The 
structure of the paper is as follows. Section [TT1 presents 
the background and the linear perturbations of the in- 
teracting dark matter-dark energy model explored here. 
Sections IIIII and IIVI describe the CMB lensing extraction 
method and the future CMB data simulation used in our 
numerical analysis, respectively. We present our results 
in Sec. W\ and draw our conclusions in Sec. IVI1 



II. PRELIMINARIES 

At the level of the background evolution equations, 
one can generally introduce a coupling between the dark 
matter and dark energy sectors as follows: 



p dm + ST-Lpdm = aQ , 
pde + ST-Lpdei^ + w) = -aQ, 



(1) 
(2) 



In this paper we investigate how allowing for a feasi- 
ble interacting dark matter and dark energy model will 



where the bars denotes background quantities, p~dm(pde) 
refers to the dark matter (dark energy) energy density, 
the dot indicates derivative with respect to conformal 
time dr = dt/a and w = Pde/pde is the dark-energy 
equation of state (P denotes the pressure). We take 
% = a/ a as the background expansion rate. We work 
in the context of a Friedman-Robertson- Walker (FRW) 
metric, assuming a flat universe and pressureless dark 



2 



matter Wd m = Pdm/pdm — 0. Q encodes the dark cou- 
pling and drives the background energy exchange be- 
tween dark matter and dark energy. In order to deduce 
the evolution of the background as well as the density 
and velocity perturbations in coupled models, we need 
an expression for the energy transfer at the level of of 
the stress-energy tensor: 

V„T$ = Ql a) , (3) 

where a = dm, cfe, refers to the energy-momentum 
tensor and is the energy-momentum transfer be- 

tween the dark matter and dark energy fluids. We con- 
sider 

Q(*"> = tU P deut n) = -Qi de) , (4) 

where £ is a dimensionless coupling (considered constant, 
as well as w, in the present analysis). H and pa e re- 
fer to the total expansion rate and dark energy density, 
background plus perturbation, i.e H — %/a + 5H and 
Pde — Pde + &Pde respectively. In the previous Eqs. ([T]) 
and Q corresponds to ^Hpde/a, see our choice of cou- 
pling in Eq. (@}. Notice from Eq. ^ that has been 
chosen parallel to the dark matter four velocity u^ m \ 
in order to avoid momentum transfer in the rest frame 
of the dark matter component [16]. For this choice of 

energy exchange Qi a \ positive (negative) values of the 
coupling £ will lead to lower (higher) dark matter energy 
densities in the past than in the uncoupled £ = case. 
In the following, we restrict ourselves here to negative 
couplings and w > 1, which avoids instability problems 
in the dark energy perturbation equations, see Ref. [l9j ] . 

The interacting model given by Eq. Q has already 
been previously explored under several assumptions, see 
Refs. fl6l - fl9l |. In those works, the linear perturbation 
analysis did not include perturbation of the expansion 
rate 5H. Let us mention that the former is included in 
the numerical analysis presented here. The latter is quite 
relevant for the correct treatment of gauge invariant per- 
turbation but it does not affect much the physical results. 
Details of the complete linear perturbation analysis will 
be presented in Ref. [24j including the specification of 
the initial conditions which have been chosen adiabatic 
for all the components except for the dark energy fluid, 
see Ref. (24|. For the numerical analysis presented here, 
we have modified the publicly available C AMB code [26[ , 
taking into account the presence of the dark coupling in 
both the background and the linear perturbation equa- 
tions. 



III. LENSING EXTRACTION 

The analysis presented here includes, in addition to 
the primary CMB anisotropy angular power spectrum, 
the information from CMB lensing. Gravitational CMB 
lensing, as already shown (see e.g. [23l [25j. can improve 



significantly the CMB constraints on several cosmolog- 
ical parameters, since it is strongly connected with the 
growth of perturbations and gravitational potentials at 
redshifts z < 1 and therefore, it can break important de- 
generacies. The lensing deflection field d can be related 
to the lensing potential <f> as d = V<p [27[. In harmonic 
space, the deflection and lensing potential multipoles fol- 
lows 

of = -tVia+Tfap, (5) 

and therefore, the power spectra Cf d = (d^df 1 *) and 
Of* = (0f are related through 

Cf d = 1(1 + l)Cf . (6) 

Lensing introduces a correlation between different 
CMB multipoles (that otherwise would be fully uncor- 
related) through the relation 

(aTbf) - (-l) m SZ'Sl'ct + SIT'/Vf , (7) 

LM 

where a and b are the T, E, B modes and 5 is a linear 
combination of the unlensed power spectra Cf h (see [28[ 
for details). 

In order to obtain the deflection power spectrum from 
the observed C ; af> , we have to invert Eq. ([7]), defining a 
quadratic estimator for the deflection field given by 

d(a,b)f!=nf W( ai b)T?i?a?W\ (8) 

W mm' 

where n a L h is a normalization factor needed to construct 
an unbiased estimator (d(a, b) must satisfy Eq. ©). This 
estimator has a variance: 

(rf(a ) 6)f*rf(a',6')^'> = <5f<(Cf + N f^') ( 9 ) 

that depends on the choice of the weighting factor W 
and leads to a noise 7V£ a bb on the deflection power 
spectrum Cjf obtained through this method. In the 
next section we describe the method followed to extract 
the lensing noise. 



IV. FUTURE CMB DATA ANALYSIS 

We evaluate the achievable constraints on the coupling 
parameter £ by a COSMOMC analysis of future mock 
CMB datasets. The analysis method we adopt here is 
based on the publicly available Markov Chain Monte 
Carlo package cosmomc 29] with a convergence diagnos- 
tic using the Gelman and Rubin statistics. We sample 
the following seven-dimensional set of cosmological 
parameters, adopting flat priors on them: the baryon 
and cold dark matter densities ili,h 2 and f2 c /i 2 , the ratio 
of the sound horizon to the angular diameter distance at 
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decoupling 9 S , the scalar spectral index n s , the overall 
normalization of the spectrum A s at k = 0.002 Mpc -1 , 
the optical depth to reionization r, and, finally, the 
coupling parameter £. 

We create full mock CMB datasets (temperature, 
E-polarization mode and lensing deflection field) with 
noise properties consistent with Planck [30j and EPIC [? 

experiments, see Tab. U for their specifications. 
The fiducial model is chosen to be the best-fit from 
the WMAP analysis of Ref. [J with fl b h 2 = 0.0227, 
n c h 2 = 0.113, n s = 0.963, r = 0.09 and £ = 0, fixing 
w = —0.9 for our numerical calculations. 



deflection field Cf d 

Kfdd \ 

I ST (\raba'b'\-l 1 V > 

l^aa'bb' l iV I I 

We compute the minimum variance lensing noise for both 
Planck and EPIC experiments by means of a publicly 
available routine [Hj]. 

The datasets (which include the lensing deflection 
power spectrum) are analyzed with a full-sky exact like- 
lihood routine 1321. 



V. RESULTS 



Experiment Channel FWHM AT/T 



Planck 


70 


14' 


4.7 




100 


10' 


2.5 




143 


7.1' 


2.2 


fsky = 0.85 








EPIC 


70 


12' 


0.05 




100 


8.4' 


0.05 




150 


5.6' 


0.06 


fsky = 0.85 









TABLE I: Planck and EPIC experimental specifications. 
Channel frequency is given in GHz, FWHM (Full- Width at 
Half-Maximum) in arc-minutes, and the temperature sensi- 
tivity per pixel in jj,K/K. The polarization sensitivity is 
AE/E = AB/B = V2AT/T. 



Parameter 


Planck 


Planck Lens 


EPIC 


EPIC Lens 


A(Q b h 2 ) 


0.00015 


0.00012 


0.00004 


0.00003 


A(n c h 2 ) 


0.0297 


0.0296 


0.0295 


0.016 


A(6 S ) 


0.00229 


0.00215 


0.00216 


0.00102 


A(r) 


0.0047 


0.0041 


0.0022 


0.0021 


A(n B ) 


0.0037 


0.0029 


0.0018 


0.0013 


A(log[10 10 A s ]) 


0.016 


0.012 


0.005 


0.004 


A(H ) 


2.34 


2.16 


2.17 


1.33 


A(fi A ) 


0.064 


0.061 


0.062 


0.038 


t 


> -0.59 


> -0.54 


> -0.56 


> -0.34 



TABLE II: 68% c.l. errors on cosmological parameters. Upper 
limits on £ are 95% c.l. constraints. 



We consider for each channel a detector noise of 
ur 1 = (9a) 2 , where 6 is the FWHM (Full- Width at Half- 
Maximum) of the beam assuming a Gaussian profile and 
a is the temperature sensitivity AT/T (see Tab.Q]for the 
polarization sensitivity). We therefore add to each Cg 
fiducial spectra a noise spectrum given by: 

N e = w- x exp(l{l + l)/l%), (10) 

where l b is given by l b = y/& In 2/9. 

In this work, we use the method presented in [28j to 
construct the weighting factor W of Eq. ([8]). In that pa- 
per, the authors choose W to be a function of the power 
spectra Cf h , which include both CMB lensing and pri- 
mary anisotropy contributions. This choice leads to five 
quadratic estimators, with ab = TT, TE, EE, EB, TB; 
the BB case is excluded because the method of Ref. [28[ 
is only valid when the lensing contribution is negligible 
compared to the primary anisotropy, assumption that 
fails for the B modes in the case of Planck. In the case of 
EPIC we have decided to neglect the BB channel since 
it may be contaminated by unknown foregrounds and/or 
it may be used for foregrounds removal. The results pre- 
sented here for the EPIC mission have therefore to be 
considered as conservative. The five quadratic estima- 
tors can be combined into a minimum variance estimator 
which provides the noise on the power spectrum of the 



Parameter 


Planck Lens 


CMBPOL Lens 


A(n b h 2 ) 


0.00013 


0.00003 


A(n c h 2 ) 


0.0010 


0.0003 


A(0a) 


0.00026 


0.00005 


A(r) 


0.0042 


0.0022 


A(n s ) 


0.0031 


0.0014 


A(/o ff [10 10 J 4s]) 


0.013 


0.005 


A(H ) 


0.53 


0.12 


A(n A ) 


0.005 


0.001 



TABLE III: 68% c.l. errors on cosmological parameters from 
Planck and EPIC with lensing extraction in the standard non 
interacting case (£ = 0). 

Table [TTI summarizes the errors from Planck and EPIC 
future data on the main cosmological parameters when a 
coupling £ among the dark energy and dark matter fluids 
is introduced in the model. Notice that the errors on the 
cosmological parameters which are degenerate with £ are 
larger than the errors that one would get for these pa- 
rameters within a standard cosmology, where £ = (see 
Tab. IIII[) . In particular, in the case of Q c h 2 , Hq, tt\ and 
9 S , the errors we obtain here are one order of magnitude 
larger than the ones obtained with £ = 0. For com- 
parison, we show in Tab. |H] the parameter constraints 
both with and without lensing extraction. The cosmo- 
logical parameter constraints from EPIC mock data are 
stronger than those coming from future Planck data when 
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the CMB lensing signal is exploited. The reason for that 
is because the EPIC experiment is expected to drastically 
reduce the noise in the CMB lensing extraction. 

Temperature power spectra 



u 

- 



+ 




2500 3000 



FIG. 1: Temperature power spectrum signal plus noise for 
Planck (top panel) and EPIC (bottom panel) experiments. 
The black curve depicts the ACDM model with Q. c h 2 = 0.113. 
The red curve illustrates a coupled model allowed by Planck 
data, with £ = -0.4 and Sl c h 2 = 0.0463. 



Lensing deflection power spectra 



- ACDM with Planck noise 

- %—OA £2 1^=0.04630 




FIG. 2: Lensing deflection power spectrum signal plus noise 
for Planck (top panel) and EPIC (bottom panel) experiments. 
The black curve depicts the ACDM model with Q. c h 2 = 0.113. 
The red curve illustrates a coupled model allowed by Planck 



data, with £ = —0.4 and Q c h 



0.0463. 



Figures (JTJ and (J3J illustrate the CMB temperature 
and lensing deflection spectra plus noise from Planck and 
EPIC experiments assuming two possible cosmologies: a 
ACDM universe with fl c h 2 = 0.113 and a coupled model 
with £ = —0.4 and fl c h 2 = 0.0463. From measurements 
of the CMB unlensed temperature spectrum (see Fig. (JTJ) 
, these two models will be degenerate, since they have 
identical spectra, albeit with different cold dark matter 



densities, since a more negative coupling can be com- 
pensated with a lower Q c h 2 . As already pointed out in 
Ref. [l9j], in a universe with a negative coupling £, the 
matter content in the past can be higher than in the 
standard ACDM scenario due to an extra contribution 
proportional to the dark energy component, and there- 
fore Q, c h 2 is strongly correlated with the coupling £. 

Notice from Fig. (JTJ that neither Planck nor EPIC data 
will be able to distinguish among coupled and uncou- 
pled models using only primary CMB anisotropy data. 
However, these two models predict distinguishable lens- 
ing potential spectra, see Fig. (J2J, and, while the Planck 
experiment will not have enough sensitivity to distinguish 
coupled versus uncoupled models, the EPIC experiment, 
with a greatly reduced noise on the CMB lensing extrac- 
tion, will be able to test coupled models. 




0.75 0.8 0.85 0.9 0.95 -0.6 -0.4 -0.2 



0.02 0.04 0.06 0.08 0.1 



FIG. 3: The panels show the 68% and 95% confidence 
level contours combining the five most correlated parameters 
(n c h 2 , 9 S , H , SI A and £) arising from a fit to mock Planck 
data without including lensing extraction in the analysis. 




FIG. 4: The panels show the 68% and 95% confidence 
level contours combining the five most correlated parameters 
{Q c h 2 , 9 S , H , fl A and f) arising from a fit to mock EPIC 
data without including lensing extraction in the analysis. 
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0.02 0.04 0.06 0.08 0.1 0.02 0.04 0.06 0.08 0.1 1.052 1.054 1.056 1.058 

n h n h e 




0.75 0.8 0.85 0.9 0.95 -0.6 -0.4 -0.2 0.02 0.04 0.06 0.08 0.1 

Q A I "c h 



FIG. 5: The panels show the 68% and 95% confidence 
level contours combining the five most correlated parameters 
(fi c /i 2 , 6> s , Ho, n A and £) arising from a fit to mock Planck 
data including lensing extraction in the analysis. 




0.06 0.08 0.1 1.052 1.054 1.056 -0.4 -0 2 

n ii a ^ 




0.75 0.8 0.85 -0.4 -0.2 0.06 0.08, 0.1 



FIG. 6: The panels show the 68% and 95% confidence 
level contours combining the five most correlated parameters 
{n c h 2 , e s , H , Q A and £) arising from a fit to mock EPIC 
data including lensing extraction in the analysis. 



Figures © and dU depict the 68% and 95% confidence 
level contours combining the five most correlated param- 
eters arising from a fit to mock Planck and EPIC data re- 
spectively, without considering CMB lensing extraction. 
We can see that the cold dark matter density £l c h 2 , the 
Hubble constant Hq, the cosmological constant the 
sound horizon angle 9 S and the coupling £ are all strongly 
correlated and neither Planck nor EPIC data will be able 
to break these degeneracies. Notice as well that, despite 
the technological advances of EPIC, the error on the cos- 
mological parameters achieved by the Planck experiment 
will not be further improved by EPIC data if no lensing 
signal is considered. 

Figures (|5|) and show the 68% and 95% confidence 
level contours combining the five most correlated param- 
eters arising from a fit to Planck and EPIC future data 
respectively, considering the information from CMB lens- 



ing. Notice that the inclusion of lensing power spectrum 
improves drastically the EPIC constraints. However, the 
addition of CMB lensing information does not change 
Planck results. This is due to the fact that the lensing 
noise for EPIC is significantly lower than for the Planck 
experiment, and therefore EPIC data would be able to re- 
ject models that otherwise would be accepted by Planck, 
see Fig. ©. 



VI. CONCLUSIONS 

The current accelerated expansion of the universe is 
driven by the so-called dark energy. This negative pres- 
sure component could be interpreted as the vacuum en- 
ergy density, or as a cosmic, dynamical scalar field. If 
a cosmic quintessence field is present in nature, it may 
couple to the other fields in nature. While the couplings 
of the quintessence field to ordinary matter are severely 
constrained, an energy exchange among the dark matter 
and dark energy sectors is allowed by current observa- 
tions. 

The major goals of the on-going Planck and the future 
EPIC experiments are to determine the nature of the 
dark energy component and to measure the remaining 
cosmological parameters with unprecedented precision. 
Several studies in the literature have been devoted to ex- 
plore the performance of Planck and EPIC experiments 
in the dark energy scenario, see for instance Ref. [23j . In 
this paper we have explored the performance of Planck 
and EPIC experiments in alternative dark energy cos- 
mologies, more concretely, in a universe with a coupling £ 
among the dark energy and dark mater components [19j . 

We have generated mock data for the Planck and EPIC 
experiments. CMB gravitational lensing extraction has 
also been included in the analysis. The lensing noise has 
been computed by means of the minimum variance es- 
timator method of Ref. [28]. The mock data have then 
been analyzed using MCMC techniques to compute the 
errors on the several cosmological parameters considered 
here. We find that relevant degeneracies are present 
among the coupling £ and some other cosmological pa- 
rameters, as the cold dark matter density ft c h 2 . There- 
fore, in the presence of a coupling, the expected Planck or 
EPIC errors on quantities as the cold dark matter energy 
density or the angular diameter distance at decoupling 8 S 
are one order of magnitude larger than in standard cos- 
mologies with £ = 0. 

When gravitational CMB lensing extraction is included 
in the analysis, Planck results remain unchanged, due 
to the high level of lensing noise for this experiment. 
However, the EPIC mission, which will benefit from a 
much lower lensing noise level, can (a) provide tighter 
constraints on the cosmological parameters, even in the 
presence of a coupling, and (b) distinguish among cou- 
pled and uncoupled models that would look identical if 
they were fitted to Planck (lensed or unlensed) data. 
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